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Abstract 

The spin-1 quantum Ising systems with three-spin interactions on two-dimensional triangular 
lattices are studied by mean-field method. The thermal variations of order parameters and phase 
diagrams are investigated in detail. The stable, metastable and unstable branches of the order 
parameters are obtained. According to the stable conditions at critical point, we find that the 
systems exhibit tricritical points. With crystal field and biquadratic interactions, the system has 
rich phase diagrams with single reentrant or double reentrant phase transitions for appropriate 
ranges of the both parameters. 
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I. INTRODUCTION 



The models with multispin interactions have attracted much attentions, such as the eight- 



vertex model 



2[ which involves both two- and four-spin interactions, the Ising model with 



four-spin interactions 



and the Baxter- Wu model js] which is a non-trivial Ising model 
with three-spin interactions on a triangular lattice. The Hamiltonian of the Baxter- Wu 
model is 

H = -J^aiajak, (1) 

(ijk) 

where (ijk) denotes a triplet of nearest neighbor sites of a triangular lattice, J is the cou- 
pling constant, and (Xj = ±1 is the Ising spin located at site i. Besides the two-dimensional 
Ising model [6], the Baxter- Wu model is another exactly solvable one which exhibits sponta- 
neous symmetry breaking. The system has been studied by series-expansion analysis [t-O], 
Monte Carlo Monte Carlo renormalization group 12, 13| and conformal invariance ap- 

proaches 15|, etc. It has the same critical tem per ature as the Ising model with three-spin 



interactions in one direction on a square lattice [16|, and belongs to the same univers ality 
class as two-dimensional 4-state Potts system with the equal critical exponents 17|, Il8| . 
namely, their specific heat critical exponent and the correlation length exponent are equal 
to 2/3 {a = u = 2/3), and correlation function exponent is r/ = 1/4. In addition, Costa et 



al. 



19| have studied the spin-1 Baxter- Wu model with a crystal field and found that the 



phase diagram presents a tri critical point. 

In 1990s, due to describing incommensurate phases of quantum one-dimensional magnet- 



20| with Hamiltonian 



ics, Tsvelik constructed a modified XXZ model 

N N 

H = ^Cr„.(T„+i + J2^Crn~l ■ (0"n X 0-n+l) , (2) 
n=l n=l 

where a" (a = x, y, z) are pauli matrices, is the total number of spins, J2 is the three- 
spin interactions strength. With competing of the two and three-spin interactions, the model 



exhibits several quantum phase |2l|,|22|. In this paper, imitating Eq. ([2]) we define a quantum 



Ising model with three-spin interactions on a triangu 



It is well-known that the Blume-Capel (BC) model 23l-l25| and the Blume-Emery-Griffiths 



ar 



attice. 



(BEG) model 26| present much interesting critical phenomenon: tricritical point, multicrit 



ical phase diagrams 27H3l|. the reentrant phenomenon 31 



metastable and unstable 



2 



states 



34 



35|. Similar to the BC model and BEG model, we study the spin-1 quantum 



Ising systems with a crystal field and biquadratic interactions by mean-field method. Our 
emphases are the thermodynamic properties and phase diagrams of the systems. We expect 
our system can exhibit that interesting proprieties. 

The outline of this paper is as follows. In section HIl using the mean-field method, we 
study the spin-1 quantum Ising system with a crystal field. In Sec. Illlt the spin-1 quantum 
Ising system with biquadratic interactions is studied in detail, and Sec. lIVl is a brief summary. 



II. QUANTUM ISING SYSTEM WITH A CRYSTAL FIELD 



The spin-1 quantum Ising model with three-spin interactions in presence of a crystal field 
on a triangular lattice is defined by the Hamiltonian 

\2 



(3) 



{ijk) 



where Si = Sfi + Sfj + S^k, S'f (a = x, y, z) are components of the spin-1 operator at site 
z, and D plays the role of a crystal field. S*" (a = x, y, z) are 



S'' = — 



'^O 1 0^ 



1 1 
1 



-i 

1 -i 
i 



\ 



1 ^ 




0-1 



In the following discussion, the coupling constant J > 0, that is, a ferromagnetic case is con- 
sidered. Based on the mean-field approximation {sg], |33], the three-spin cluster Hamiltonian 
iJi23 can be written as 



H 



123 



J SI ■ (S2 X S3 ) +D [{Slf + (^2^)' + {Slf) - 4Jm {S^ + S| + 



(4) 



where m is the average magnetization of the cluster, {i.e., m = {{Sf + 5*1 + 5*1) /3)). In the 
representation of the direct product of S^, S2 and S^, H123 can be written as a form of 
27 X 27 matrix. Thus the partition function Z = Tr[exp (— /3i^i23)] ( (3 = -j^, where is 
the Boltzmann constant, T denotes the temperature of the system) can be obtained as 

Z = 1 + 4ai cosh{4:Kmf + 4a2 cosh (AKm) + 2e'^^^" cosh {12Km) , (5) 

in which 

= e-2^^" (1 + 2 cosh ir), (6) 



a2 = e-2^^« 2 cosh (KDq) + cosh \^K^J2 + D^j , (7) 

where K = J/ksT is the reciprocal of the reduced temperature, and Dq = D/J is the 
reduced crystal-field parameter. The reduced free energy / of the three-spin cluster in the 
system as a function of K and m is further given by 

/ (m, K) = Grri^ - K'^ In [Z (m, K)] . (8) 

According to the equilibrium condition of the system {df /dm = 0), we can obtain the 
self-consistent equation of the magnetization 

m = — , (9) 

3Z' ^ ^ 

where 

a = 2ai sinh {8Km) + 2a2 sinh (AKm) + 3e~^^^° sinh {12Km) . (10) 

By solving the self-consistent Eq. (Q, numerically, we obtain the thermal variations of 
the magnetization for several values of Dq (Fig. 1). We obtain the stable, metastable and 
unstable branches, which are represented by the solid, dot and dash-dotted lines, respectively. 
It can be seen that for Dq = 1, 1.6 and 1.9 the stable magnetization are continuous, indicating 
the second-order phase transitions. For Dq = 1.6 and 1.9, we also find unstable solutions. On 
the other hand, the stable magnetization drops discontinuously from a finite value nip (Dq) 
to zero at a temperature Tp {Dq), which characterizes a first-order phase transition. The 
first-order phase transition temperature Tp (-Do) can be obtained by the free energy, when 
its local minimum at m 7^ is equal to the local minimum at m = 0, Is^]- When the 
reduced crystal-field parameter is slightly larger than the tricritical point D}^ (1.95138), the 
system can undergo the first-order phase transition, such as, for Dq = 1.98 the first-order 
phase transition temperature is ksTp/J ^ 0.84112. 

Figure 2 shows the phase diagram of the spin-1 quantum Ising system in the presence of 
a crystal field. The critical line separates the ferromagnetic phase from the paramagnetic 
phase. We can see that the second-order phase transition (solid line) occurs for finite values 
of the crystal field. The second-order phase transition temperature falls smoothly with the 
reduced crystal- field parameter Dq increasing and reaches the tricritical point. The value of 
the tricritical point is Dq = 1.95138 and ksTtj-fJ = 1.1898, which is different from that in 
ref [13] D^ = 1.3089 and kBT^^/ J = 1.2225. Then, the first-order phase transition line falls 



smoothly from the tricritical point to zero. This phenomenon is similar to the spin-1 Baxter- 
Wu model [l^. It need to clear that the last point of the curve in Fig. 2 is Do = 2.00524999 
and kBTp/J = 0.4635712. 

III. QUANTUM ISING SYSTEM WITH BIQUADRATIC INTERACTIONS 
A. Model and formulation 

The Hamiltonian of the spin-1 quantum Ising system with biquadratic interactions on a 
triangular lattice is defined as 

H = -jJ2s!- {s, X s^y -Gj^istf {s]Y + dJ2 {Stf , (11) 

{ijk) (ij) i 

where (ij) indicates all nearest-neighbor pairs of sites of a triangular lattice and G denotes 
the biquadratic interactions parameter. Different from that in Sec. [Tll for solving this model, 
we need two order parameters: the magnetization m and quadrupolar moment q. Using 
the similar method as in Sec. [Til the three-spin cluster Hamiltonian -^123 has the following 
expression: 

^123 = -J St ■ [S, xS,y-G {{Srf {Slf + {S^,f {Slf + {Slf {Slf) (12) 
+ {D- AGq) {{Slf + (5|)' + [Slf) - 4Jm {SI + S| + ^3^) , 

where m has been defined in Sec. [TTl <? is the average squared magnetization of the three- 
spin cluster, q = \ [{Sl) + (S"!) + (S"!) ) /3). According to the definition of the partition 
function, it is obtained as 

Z = l + 2bi cosh(12Km) + 4^2 cosh {AKmf + Ah^ cosh {AKm) , (13) 



in which 



62 = e^(-2^o+Go+8<?Go)^i + 2coshir), (15) 



^3 = V8+(-2^o + (3 + 8g)Go)', (16) 

b, = bi + e-^(^"-4'^«o) + eKi-2Do+3Go/2+S,Go) ^^^^ {Kh/2) , (17) 
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where Go = G/J is the reduced biquadratic parameter. The reduced free energy of the 
three-spin cluster can be written as 

/ {K, Go, Do, m, q) = 6m^ + GGoQ^ - R-^ In Z {K, Gq, Dq, m, q) . (18) 

With the equihbrium conditions {df /dm = and Of /dq = 0) of the system, the self- 
consistent equations of the two order parameters are 

m = 5I (19) 

and 

(-) 

in which 

g = 3bi sinh {12Km) + 262 sinh {8Km) + 264 sinh (AKm) , (21) 
h = 3bi cosh(12/s:m) + 4^2 cosh {4Km f + 2(26i + 64 + ^5) cosh {4Km) , (22) 

where 

65 = ei^(-2Do+3Go/2+8gGo)[^ogh (K63/2) + sinh (K63/2) ((3 + 8q) Go - 2Do)/h]. (23) 
B. Results and Discussions 

By numerically calculating Eqs. f|T9|) and fl20|) . the behavior of the magnetization m and 
quadrupolar moment q as functions of temperature are obtained (see Figs. 3 and 4). We 
obtain the stable (Ml and Ql), unstable (M2 and Q2) and metastable (M3 and Q3) branches. 
As in Sec. HIl the stable branches of the order parameters undergo the first-order or second- 
order phase transition. What is interesting is that in Figs. 3(d) and 3(e), the stable branches 
of the order parameters undergo one first-order phase transition and one second-order phase- 
transition successively. Fig. 4(b) shows that the stable branches of the order parameters 
undergo two second-order phase transitions and one first-order phase transition. Otherwise, 
what should be mentioned is that the stable branch of the quadrupole order parameter does 
not undergo any phase transition when the stable branch of the magnetization is equal to 
zero [35'] (see Fig. 3(f)). Seen from Fig. 3, at T = 0, the system can be in the ground state 
with m = q = 1, m = q = 2/3 or m = q = 0, which is similar to the ground state of 
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the S = 1 Ising model on triangular lattice [38|. The actual ground state depends on the 



competition between the reduced biquadratic interaction and crystal field. 

Figure 5 shows three phase diagrams for different biquadratic parameter Gq. It can 
be seen that the first-order phase transition can occur either in the range of Dq < Dq 
(Fig. 5(a)) or Dq > (Fig. 5(c)). At fixed reduced biquadratic and crystal- field parameter 
in Fig. 5(a), when the temperature decrease the transition from paramagnetic phase (P) 
to ferromagnetic phase (F) is encountered, and if the temperature is lowered further, the 
reentrant phase transition from ferromagnetic phase to paramagnetic phase can take place. 
Fig. 5(c) shows a double reentrant phenomena, that is, at fixed reduced biquadratic and 
crystal field parameter, as the temperature decrease, the P— J-F— J-F sequence of phase is 
encountered. It may be determined by the competition among quantum effect, crystal- field 
and the biquadratic interaction. 



IV. SUMMARY 



In this paper, we have investigated the spin-1 quantum Ising systems with three-spin in- 
teractions on two-dimensional triangular lattices by mean-field approximation. The thermal 
variations of the magnetization m and quadrupolar moment q are obtained. By comparing 
the free energy of these solutions, we obtain the stable, unstable and metastable branches of 
the order parameters. In addition, we find that the quantum Ising systems present tricriti- 
cal points for finite values of the crystal field, and the system with biquadratic interactions 
exhibits single reentrant and double reentrant phase transitions. 
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Fig. 1 The magnetization m versus temperature ksT/J for the spin-1 quantum Ising 
system with a crystal field on a triangular lattice. Curves a, b, c, d represent the stable 
branches for Dq — 1, 1.6, 1.9, 1.98, respectively. Curves b', c', d' represent the unstable 
branches and curve d" is the metastable branch. 

Fig. 2 The phase diagram of spin-1 quantum Ising system with a crystal field on 
a triangular lattice. The solid line represents the second-order phase transition and 
dot line represents the first-order phase transition. The big real point (solid dot) is the 
tricritical point. P and F represent the paramagnetic and ferromagnetic phases, respectively. 

Fig. 3 The magnetization m and quadrupolar moment q as functions of temperature 
ksT/J and crystal-field parameter Do with Go = —0.5 for the spin-1 quantum Ising model 
with biquadratic interactions on a triangular lattice. Ml, Ql (sohd line) are the stable 
branches of the magnetization m and quadrupolar moment q, respectively; M2, Q2 (dot 
line) represent the unstable branches; M3, Q3 (dash-dot line) represent the metastable 
branches. Tp and Tq indicate the first-order and second-order phase transition temperature, 
respectively. 

Fig. 4 The nonzero order parameters m and q as functions of temperature ksT/J and 
crystal-field parameter Dq with Go — S for the spin-1 quantum Ising model with biquadratic 
interactions on a triangular lattice. Ml, Ql (solid fine) represent the stable branches; 
M2, Q2 (dot fine) represent the unstable branches. Tp and Tc indicate the first-order 
and second-order phase transition temperatures, respectively. is the unstable critical 
temperature. 

Fig. 5 Phase diagrams of the spin-1 quantum Ising model with biquadratic interactions 
on a triangular lattice for different Go- The solid and dot lines represent the second-order 
and first-order phase transition lines and the solid dot represents the tricritical point. P and 
F denote the paramagnetic and ferromagnetic phase, respectively, (a) shows a P— >F— >P 
reentrant phase transitions, (c) shows a P— >F— >P— >F double reentrant phase transitions. 
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